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Abstract. We provide an informal overview on the theory of transport equations with 
non smooth velocity fields, and on some applications of this theory to the well-poscdness of 
hyperbolic systems of conservation laws. 



1. Introduction 

Our attempt in this note is to give an overview on the theory of transport and continuity 
equations 

d t u + b-Vu = and d t u + div (bu) = , where b : [0, +oo[xR d -> R d , 

and on its applications to some classes of nonlinear PDEs. Since various extensive surveys 
and lecture notes on this topic are available (for instance the interested reader could consult 
0, IH dU [T3], where also many references to the literature are provided), we prefer here to 
follow a more informal approach, by starting from the applications and then by step by step 
addressing the related issues concerning transport and continuity equations. 

We now describe the topics that shall be discussed in this note: the related references will 
be provided in the corresponding sections. 

We start in Section [2] by going over some of the main aspects of the theory of conservation 
laws. For this class of PDEs, which typically model conservation of physical quantities, ex- 
istence fails in the class of classical (smooth) solutions, but bounded distributional solutions 
are not unique. Admissibility conditions have been designed in the attempt of restoring 
uniqueness of distributional solutions: an example is the entropy dissipation condition (j2.3|) . 
namely the requirement that some nonlinear function of the solution must be dissipated. 
The underling idea is that the entropy criterion should select the physically relevant solu- 
tion. General existence and uniqueness results concerning conservation laws are available in 
the scalar case, and for one-dimensional systems under smallness assumptions on the total 
variation of the initial data. However, proving existence and uniqueness results valid for 
general systems in several space dimensions is still a major open problem. 

We then turn our attention to one particular multi-dimensional system, the so-called 
Keyfitz and Kranzer system, see (12.51) . The peculiarity of (12.51) lies in the dependence of the 
nonlinearity on the modulus of the solution only, thus inducing a sort of radial symmetry. 
Formally, this system splits into the scalar conservation law (12.61) and the transport equation 
(12. 7p . and this paves the way for the application of transport equation techniques to the 
study of (T2~5|) . 
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However, the only regularity (with respect to the space variable) we can expect on the 
velocity field appearing in (12 .7p is BV regularity. Hence, we are forced to discuss the be- 
haviour of this transport equation in the case of non smooth velocity fields. We go over this 
research area in Section [31 addressing some of the main points of the well-posedness theory 
in the cases the velocity field enjoys either Sobolev or BV regularity. 

We finally motivate in Section H] how these techniques lead to rigorous well-posedness 
results for the Keyfitz and Kranzer system (12.51) : actually, one has to introduce some ad- 
ditional arguments which take care of the fact that no information on the divergence, but 
only on the compression rate of the velocity fields are typically expected in this context. 
We conclude the note by collecting some results, analogue in spirit but valid in one space 
dimension, regarding again the Keyfitz and Kranzer system, and also the chromatography 
system ( 14.61) . 

2. Systems of conservation laws 

In this section we informally go over some results concerning systems of conservation laws 
in the form 

d 

d t U(t,x)+J2d*MU(t,x)] =0. (2.1) 
i=i 

The unknown U is a function 

U(t,x) : [0,+oo[xM d ^R fe 

and all the functions Fi, i — 1, . . . , d are defined on an open set of lR fc and take values in IR fc . 
Classical references concerning these topics are the books by Serre [30] and Dafermos |15j . 
Also, the book by Bressan [9] focuses on the study of the Cauchy problem for systems in one 
space dimension, d — 1, while the recent book by Benzoni-Gavage and Serre [8] is concerned 
with systems in several space dimensions. 

Conservation laws like (12. ip have several and diverse physical applications. In particular, 
the Euler equation can be written in the form (12. ip . The Euler equation is obtained from 
the Navier-Stokes equation by formally setting the viscosity coefficients equal to zero and 
describes the motion of the particles in an inviscid fluid. 

One of the main difficulties in the study of (12.11) is the possible loss of regularity of the 
solution. More precisely, let us first focus on the Cauchy problem obtained by coupling (12.11) 
with the initial condition 

U(0,x) = U(x). (2.2) 

As a matter of fact, even if the datum (12.21) is very regular, in general there is no classical 
solution of (12.11) defined on the whole time interval t G [0, +oo[. Here a classical solution is 
a function U G C 1 ([0, +oo[xIR ,i ; M fc ) satisfying (I2.1l) - (l2.2p pointwise. Remarkably, examples 
of solutions starting from a datum U G C°° and developing discontinuities in finite time are 
available even in the scalar, one- dimensional case k = d = 1. 

In the attempt of obtaining global existence results, it is thus natural to introduce a 
notion of weak solution, interpreting (12.11) in the sense of distributions. However, in general 
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a weak solution of the Cauchy problem (I2.ip - (I2.2I) is not unique: as for the development of 
discontinuities, examples of non uniqueness for a distributional solution of (12.11) are available 
even in the scalar, one-dimensional case k = d = 1. In the attempt at selecting a unique 
solution, several admissibility conditions can be introduced (entropy inequalities, Lax and 
Liu conditions on shocks), often motivated by physical considerations: we refer to the book 
by Dafermos [T5l Chapters 4 and 8] for a detailed exposition. 

Since we will need it in the following, we recall the definition of entropy here. Let 
7] : R h -»• R, Q : R k -»• R d be two smooth functions, Q = (Q 1} . . .,Q d ). Then (rj,Q) is 
an entropy-entropy flux pair for (12.11) if 

V?7 • DFi = VQi for every % — 1, . . . , d. 

Here DFi denotes the Jacobian matrix of Fi and V77 the gradient of r\. A distributional 
solution of (12.11) is entropy admissible if the inequality 



holds in the sense of distribution for every (77, Q) entropy-entropy flux pair such that 77 is 
convex. 

In the scalar case k = 1, the fundamental work by Kruzkov [25] established existence and 
uniqueness of global in time, entropy admissible solutions of the Cauchy problem (I2.ip - (I2.2I) 
under the assumption that the initial datum U is bounded, U G L°°(R). Also, Kruzkov [25] 
showed that the BV regularity is propagated, namely U(t, •) G BV(R) for every t > if 
U G BV(R). 

Let us now consider the one- dimensional case, d = 1, k > 1: in this case, (12. ip becames 



The existence of a global in time, entropy admissible solution of the Cauchy problem (12.41) - 
(12. 2p was established by Glimm [22]. Uniqueness results were obtained in a series of papers 
by Bressan and several collaborators and are discussed in Bressan's book [9J. All these 
results concerning the Cauchy problem (I2.4l) - (l2.2p hold under the assumptions that the total 
variation of the initial datum U is small enough and that the system satisfies an hypothesis 
of so-called strict hyperbolicity (see for example Dafermos [151 Sections 3.1 and 9.6] for the 
exact definition and for a related discussion). If this is not the case, then either existence 
or uniqueness of global in time, entropy admissible solutions may fail (see e.g. the examples 
discussed in Jenssen [23] and Panov |27j). 

Finally, we consider the general case of a systems of conservation laws in several space 
dimensions (12 .4p . d > l,k > 1. Presently, no general result is available concerning either 
existence or uniqueness of global in time distributional solutions. Tackling these issues in 
big generality is certainly regarded as one of the most challenging open problems in this 
research area. 

For specific classes of systems, however, well-posedess results are available. In the fol- 
lowing, we focus on cases when these results are obtained by employing transport equation 
techniques. 



d t [r)(U)]+div [Q(U)] <0 



(2.3) 



d t U + d x [F(U)] = 0, where U G R k . 



(2.4) 
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The Keyfitz and Kranzer system was introduced in [23] and takes the form 





k U{0,x) = U{x). 

Here, f a : R — > R is a smooth function for each a — 1, ... ,d. By setting p —\U\ and 

= (0 1 ,...,0 fc ) = ^/|^|, 
system ( 12.51) formally splits as the coupling between the scalar conservation law 

<V + div(/(p)p) =0 



(2.6) 



and the transport equations 

Mi + M ■ V0i = o 
e i (o,x) = u(x)/\u(x)\ 



i = l,...,k. (2.7) 



Here, / = {fi, ■ ■ ■ , fd)- Some remarks are here in order: equation (12. 6p is a scalar con- 
servation law in several space dimensions and does not depend on 9. We can thus apply 
Kruzkov's [25] theory to obtain existence and uniqueness of global in time, entropy admis- 
sible, distributional solution. Then, we can plug the function p in the second equation, 
obtaining a transport equation depending on 9 only. However, even if the initial datum 
p(0, x) = \U(x)\ is extremely regular, the best we can hope is that the BV regularity is prop- 
agated, i.e. that pit,-) £ BV{R) for every t > if p(0, ■) £ BV(M). Indeed, as mentioned 
before, a solution having very regular initial data can in general develop discontinuities in 
finite time. 

In this way, we find that the analysis of (I2.6p - (I2.7I) has close links to the study of continuity 
and transport equations with weakly differentiable coefficients: an informal overview on these 
issues is provided in Section [31 while in Section H] we come back to the Keyfitz and Kranzer 
system (12.51) and we discuss some well-posedness results obtained by pursing the approach 
we sketched before. 



3. Transport equation 
In this section we are concerned with the well-posedness of the transport equation 

d t u{t, x) + b{t, x) ■ Wu{t, x) = 
u(0, x) = u(x) 
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and of the related continuity equation 

d t u(t, x) + div (b(t, x)u{t, x)) =0 
u(0, x) = u{x) , 

in which b : [0, +oo[xIR d — > M. d is a vector field, u G L oc (M d ) is the initial datum, and the 
unknown u belongs to L°°([0, +oo[xIR d ). 

The well-poseness theory for these equations is classical and well understood in the case 
when b is sufficiently smooth (at least Lipschitz with respect to the spatial variable, uniformly 
with respect to the time), and is strongly based on the so-called theory of characteristics, 
i.e. on the connection between (13.11) and (13.21) and the ordinary differential equation 

'*«> ,3. 3) 
7(0) = X . 

However, in many applications motived by physical models, non-smooth vector fields show 
up as velocity fields of transport or continuity equations. In addition to the case of the 
Keyfitz and Kranzer system (12. 5p described at the end of Section [2J the theory of non-smooth 
transport equations has interesting applications to the study of the Boltzmann equation and 
of the Vlasov-Poisson equation: these applications were investigated by DiPerna and Lions 
in [19] and p2] respectively. 

This motivates the great interest arisen in the study of (13. ip and (13. 2p when b is only 
in some classes of weak differentiability. The first seminal result in this direction is due to 
DiPerna and Lions [20], who deal with the case of Sobolev regularity, under boundedness 
assumptions on the spatial divergence of the vector field. In [20J the notion of renormalized 
solution to the transport (or continuity) equation is introduced. Let us explain heuristically 
the motivation for such a notion. Assume for the moment that the vector field b is compactly 
supported and that the spatial divergence 

div6 G L 1 t (L^)=L 1 ([0,+oo[ t ;L oo (Ri)). (3.4) 

We multiply our equation 

d t u + b ■ Vu = (3.5) 
by 2m, and by formally applying the chain-rule we deduce 

d t (u 2 ) +b- V(w 2 ) =0. (3.6) 
We now integrate over the space M. d for every fixed valued t of the time, obtaining 



— / u 2 (t,x)dx = I div6(£, x)u 2 (t, x) dx 

< \\divb(t, •)l|x,° o (R tl ) / u 2 (t,x)dx. 



(3.7) 



thanks to the divergence theorem. Now, a simple application of Gronwall's lemma implies 
that, if the L 2 (M. d ) norm of the solution vanishes at the initial time, then it vanishes for all 
time. By the linearity of (13. ip we obtain uniqueness for the Cauchy problem. 
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In the above argument, the delicate point is hidden in the passage from (13.51) to (13. 6p . 
Indeed, since no regularity beyond boundedness is assumed on the solution u, the applica- 
tion of the chain-rule formula is not justified. This led DiPerna and Lions [20] to define 
renormalized solution of (13.11) as bounded distributional solutions u for which 

d t ((3(u))+b-V((3(u))=0 (3.8) 

holds for any function (3 G C 1 (R; R). In some sense, the validity of the "chain-rule formula" 
is assumed "by definition" for renormalized solutions. 

By relying on an argument similar to the one we gave before, one can show that unique- 
ness results are satisfied by the solutions of the transport equation (13.11) and the continuity 
equation (13. 2p . provided that the vector field b satisfies the following property: all bounded 
distributional solutions of (13.51) are renormalized and satisfy some strong continuity require- 
ment. 

The core of the argument in DiPerna and Lions [20] is thus the proof of the validity of 
(13. 8p for Sobolev vector fields, and this is achieved via a regularization procedure combined 
with a control on the convergence of the error term that appears. This gives well-posedness 
for (13.11) (and similarly for (13.21) ) assuming that b is a bounded vector field such that 

beLKW 1 /) and dwbeLl(L^). (3.9) 

More recently, this result has been extended by Ambrosio [TJ to the case of bounded vector 
fields with bounded variation with respect to the space variable (i.e. such that the space 
distributional derivative is a locally finite measure), that is one replaces (13. 9p with 

b G L\ (BV X ) and div b G L\ (Lf) . 

In both results, the need of controlling the spatial divergence of b essentially comes from 
computations analogue to (13.71) . The hard point of the proof in [1] is to control the conver- 
gence of the error term in the regularization procedure in this weaker context. The argument 
is based on more refined arguments of geometric measure theory, in particular on some fine 
properties of BV functions. For a general survey on this topic, we suggest for instance 

pi iii am . 

The possibility of extending the uniqueness result by removing the assumption b G L\(BV X ) 
is ruled out by an intriguing counterexample due to Depauw [17]: he exhibits a bounded, 
divergence free vector field such that the Cauchy problem 

J d t u + div (bu) = 
\ u(0,x) = 

admits a non trivial solution. The vector field in [UJ satisfies b G lq oc (]0; 1];BV(R 2 ; R 2 )), 
but the BV norm is not integrable up to time t — 0. 

As a side remark, we mention that a result in a recent work [5] by Ambrosio, Figalli 
and the authors establishes well-posedness for continuity equations where the vector field 
b G .BViocQO, +oo[xR d ; R d ) and has a possible blow up of the BV norm at t = 0. This does 
not contradict Depauw's counterexample because this result is obtained by imposing on the 
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solution an additional regularity assumption, namely that the map t h- > u(t, •) is continuous 
in L 1 1 oc (lR d ) endowed with the strong topology 

4. Applications to systems of conservation laws 

4.1. Nearly incompressible vector fields. Let us go back to the example of the Keyfltz 
and Kranzer system (12. 5ft : as pointed out at the end of Section [21 Kruzkov's [25] theory of 
entropy solutions ensures that the vector field 

b(t,x)=f(p)(t,x) 

appearing in the transport equation (12.71) has BV regularity with respect to the space vari- 
able. Thus, the theory described in Section [3] is concerned with the right regularity frame- 
work, i.e. with vector fields having bounded variation. However, the additional assumption 
of boundedness of the divergence of the vector field may be not satisfied, and this means 
that the theory does not apply straightforwardly to the analysis of the Keyfitz and Kranzer 
system (12. 5p . Remember that boundedness of the divergence of the vector field is apparently 
an unavoidable assumption if one tries to carry on computations like in (13.71) . 

We now illustrate some assumptions one can impose on the vector field b in (13. ip that are 
more natural in view of the applications to the Keyfitz and Kranzer system (12.51) : in this 
case, we have an additional information, coming from (12. 6p . namely that the function p is a 
BV solution of the continuity equation with vector field f(p). Also, Kruzkov's theory [25] 
guarantees that the scalar conservation law (12. 6p satisfies a comparison principle, namely 

< — < p(t, x) < C for every t > and for a.e. x G 1R 

if 1/C < p(0,x) < C for a.e. x G M. Note that requiring that p(0, x) is non negative is 
natural here, since p(0, x) is the modulus of the initial datum. 

Let us consider the (d + l)-dimensional vector field in IR/~ x defined by 

B(t, x) = (p(t, x), f(p(t, x))p(t, x)) eRxR d . 

Thus, (12. 6 p simply means that B(t,x) is divergence free in space-time, i.e. 

&iv t . x B(t,x) = . 

We add a fictitious "time variable" s G [0, +oo[, and by setting 

9i(s,t,x) = &i(t,x) 

we realize that the equation written at the first line of (12. 5p can be rewritten as 

d s i (s,t,x) + div t>x (B(t,x)0 i (s,t,x)) =0 for i = l,...,k. (4.1) 
Here we are exploiting the fact that 

e(t,x) = (ei,...,e k )(t,x) = u/\u\. 
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Thus, we can now apply the results of Section [3] to the "lifted" equation (14 .11) . since we are 
now dealing with a BV divergence free vector field (precisely thanks to the BV regularity 
of p). This means, in particular, that 8i(s,t,x) is a renormalized solution of (14.11) . that is 

d s [i3(9 i (s,t,x))]+diy ttX [B(t,x)/3(6i(s,t,x))] =0 

for all (3 G C 1 (M; M). But, &i not depending on s, this can be rewritten as 

d t [p(t,x)p(9 i (t,x))] +div [b(t,x)p(t,x)f3(e i (t,x))] =0, (4.2) 

which is reminiscent of the renormalization condition for the vector field b(t,x), apart from 
the term p(t, x) now appearing. 

The previous discussion motivates the following definition: 

Definition 4.1. We say that a bounded vector field b is nearly incompressible if there exist 
a function p and a constant C > such that 

< < p(t, x) < C < +oo for Jzf d+1 -a.e. (t, x) G ]0, +oo[xM d (4.3) 

and 

d t p + div (bp) = . (4.4) 



The renormalization property, in this context, is the requirement of the validity of (14.21) 
for every bounded function 8i solving 

t [p0j +div [bp6i] =0 

One can show that uniqueness results hold for nearly incompressible vector fields with 
the renormalization property: we refer to the notes by De Lellis [16] for a complete and 
systematic treatment of the renormalization theory for nearly incompressible vector fields 
and for various related references. The following lemma is one of the main results of this 
theory, and can be proved by relying on the "fictitious variable" argument described above. 

Lemma 4.2 ([HI Lemma 5.10]). Let 

beL°°n 5Vioc([0, +oo[xM d ; R d ) 

and 

p G L°° fl BVi oc ([0, +oo[xR d ) , with p(0, •) G BV loc (R d ) 

be such that 

d t p + div (bp) = . 
Let u G L°°([0, +oo[xR d ) be a solution of 

dt(pu) + div (bpu) = . 
Then u is renormalized, that is for every (3 G C 1 (1R; M.) there holds 

d t (p(3(u)) +div (b P p(u)) = 0. 
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A major conjecture concerns the possibility of removing the assumption of BV regularity 
on p in Lemma [4.21 

Conjecture 4.3. Any bounded nearly incompressible vector field b £ BV([0, +oo[xM d ; M. d ) 
has the renormalization property. 

This conjecture is also important in view of a compactness conjecture advanced by Bressan 
in [IT], again in connection to the question of existence of solutions to the Keyfitz and 
Kranzer system. Indeed, Ambrosio, Bouchut and De Lellis [3] proved that a positive answer 
to Conjecture 14.31 would imply a positive answer to the compactness conjecture, which reads 
as follows: 

Conjecture 4.4 (Bressan's compactness conjecture). Let : [0, +oo[xIR d — > M. d , k £ N, be a 
sequence of smooth vector fields and denote by Xk the (classical) flows associated to them. 
Assume that ||6fc||oo + ||V6fc||x,i is uniformly bounded and that the flows Xk are uniformly 
nearly incompressible, i.e. that 

— < det(V x X k (t,x)) < C for some constant C > 0. (4.5) 

Then the sequence {X k } is strongly precompact in Lj^QO, +oo[xIR d ). 

Both conjectures are presently open, apart from some particular cases of regularity. Some 
recent advances are presented in Ambrosio, De Lellis and Maly [7]. 

4.2. Keyfitz and Kranzer system in several space dimensions. By relying on Lem- 
ma [4721 and on the uniqueness results for nearly incompressible vector fields with the normal- 
ization property, Ambrosio and De Lellis [6] and Ambrosio, Bouchut and De Lellis [3] proved 
well-posedness results for the Keyfitz and Kranzer system (12.5ft . They followed the approach 
sketched at the end of Section [2J which was inspired by considerations in Bressan [TDj . 

Before stating the well-posedness result in [3], we have to introduce a preliminary defini- 
tion. Let U be a locally bounded function solving (12.51) in the sense of distributions. Then 
U is a renormalized entropy admissible solution if p — \U\ is an entropy admissible solution 
of (12.61) satisfying lim 4 ^ + p(t, •) = \U\ in the strong topology of L 1 1 oc (IR d ). By relying on a 
classification of the entropies for (12. 5p due to Panov [27J and Frid [21], under quite general 
assumptions on the function / any renormalized entropy solution is indeed an entropy ad- 
missible solution of (12.51) (for a proof of this implication see for example the notes by De 
Lellis [Ttjj Proposition 5.7]). Conversely, in general entropy admissible solutions of (12.51) are 
not renormalized and they are not unique (for an example see e.g. Bressan |lUj). 

Existence, uniqueness and stability of renormalized entropy admissible solutions are es- 
tablished in the following theorem. 

Theorem 4.5 ([3, Theorem 2.6]). Assume that the function f = (fi, ■ ■ ■ , fd) "is locally 
Lipschitz continuous and that the modulus of the initial datum satisfies \U\ £ L°°(M. d ) n 
BV\ oc (R d ). Then there exists a unique renormalized entropy admissible solution of (12.51) . 
Also, assume that U n is a sequence of initial data satisfying the following assumptions: 



10 



GIANLUCA CRIPPA AND LAURA V. SPINOLO 



(1) ||!7 n ||x,oo < C for some constant C > 0; 

(2) for every bounded open set Q C ~R d , p n = \U n \ satisfies \\p n \\BV{n) < C(fl); 

(3) U n ^U strongly in L 1 1 oc (IR o! , IR fc ). 

Then U n -> U strongly in Lj^flO, +oo[xM d , M fe ). 

As a final remark, we point out that a counterexample due to Crippa and De Lellis [H] in- 
dicates that strategies exploiting transport equation techniques aimed to show well-posedness 
for hyperbolic systems are not expected to work for general systems of conservation laws. 
More precisely, the analysis in [T4] shows that (in dimension greater than or equal to 3) there 
exists no functional space containing BV which is closed by iteration of transport equations 
and which embeds compactly in L\ oc . This suggests that approximation schemes valid for 
general conservation laws cannot be constructed by relying on transport equation techniques. 

4.3. Transport equations techniques applied to some systems in one space dimen- 
sion. If one restricts to one space dimension d — 1, well-posedness results for the continuity 
equation 

d t u + d x (bu) = 

are available under much weaker assumptions than those discussed in the previous sections. 
In particular, Panov [251 [25] proved existence and uniqueness results for bounded vector 
fields b G L°°([0, +oo[xIR) satisfying a near incompressibility condition similar to the one 
introduced in Definition 4.1. See also Serre [29]. These results are then applied to the 
analysis of the one-dimensional Keyfitz and Kranzer system 

d t U + d x [f(\U\)U] =0 

obtaining existence and uniqueness of renormalized, entropy admissible solutions with merely 
bounded initial data (the presence of linear source terms is also contemplated in [28]). See 
also Panov [27], where the vanishing viscosity approximation 



d t U E + d x [f(\U E \)U s ] =eU: 



e 

xxi 



is investigated, and the convergence of the family U e as e — > + is established. 

As an other example of applications of transport equation techniques to a system of 
conservation laws in one dimension we mention the so-called chromatography system 

, ) =° 

+ u±-\ h u k J 

(4.6) 

d t Uk + d x [ — ) = . 

+ Ml H h u k J 
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Here the unknowns Ui,...,Uk are real-valued and one is usually interesting in finding a 
solution satisfying U\ > 0, . . . , Uk > 0. By setting 

v — u\ + h Uk and Wi = Ui for i — 2, . . . , k 

one obtains that the system writes as the coupling between the scalar conservation law 

and the continuity equations 

dtWi + d x (rjA =0 for % = 2, . . . , k. (4.7) 

As in (12. 7p . the vector field b(t, x) = 1/(1 + v) in (14.71) depends on v and hence one expects 
low regularity. By following this strategy, well-posedness results for the chromatography 
system are obtained in [5]. See also the previous works by Panov [27J Remark 4 page 140] 
and Bressan and Shen [12] for similar approaches. 
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